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Abstract 



We study the low-temperature effective potential of the Ising model. We 
evaluate the three-point and four-point zero-momentum renormalized cou- 
pling constants that parametrize the expansion of the effective potential near 
the coexistence curve. These results are obtained by a constrained analysis 
of the e-expansion that uses accurate estimates for the two-dimensional Ising 
model. 
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I. INTRODUCTION 



In statistical physics the effective potential corresponds to the free-energy density T as 
a function of the order parameter, which, for spin models, is the magnetization M. The 
global minimum of the effective potential determines the value of the order parameter, 
that characterizes the phase of the model. In the high-temperature (HT) or symmetric 
phase the minimum is unique with M — 0. According to the Ginzburg-Landau theory, as 
the temperature decreases below the critical value, the effective potential takes a double- 
well shape. The order parameter does not vanish anymore and the system is in the low- 
temperature (LT) or broken phase. Actually in the broken phase the double-well shape is 
not correct because the effective potential must always be convex. In this phase it should 
present a flat region around the origin. For a discussion see e.g. Refs. [|l]|| and references 
therein. 

Recently many works have been devoted to the study of the effective potential — or, 
equivalently, of the equation of state — in the limit of small magnetization for O(N) models 
in the HT phase. Indeed, in this case, the effective potential admits a regular expansion 
around M = with coefficients that are related to the zero-momentum n-point renormal- 
ized coupling constants. Rather accurate estimates of these quantities have been obtained 
exploiting various field-theoretic approaches (see e.g. and lattice techniques (see e.g. 

mi)- 

In the LT phase the effective potential shows a complex behaviour due to the Goldstone 
bosons and does not admit a regular expansion. The only exception is the Ising model 
where the symmetry is discrete and therefore the Goldstone bosons are absent. In this case, 
if M = lim#_ > o+ M(H), for M > M (i.e. for H > 0), the effective potential admits a regular 
expansion in powers of (M — M ). The coefficients of this expansion are related to the LT 
zero-momentum n-point renormalized coupling constants g n . These quantities are harder to 
estimate than the corresponding HT couplings and indeed they are known with much less 
precision. At present there are some estimates for g^ and g±. However the uncertainty is 
often very large or difficult to estimate (especially for g 4 ) PH.H-IBI. Surprisingly one finds 



the same situation for the two-dimensional Ising model (see Refs. |]10|,|15| ). 

In this paper we study the LT effective potential for the Ising model in two and three 
dimensions. We focus on the expansion in powers of (M — Mo) that we parametrize in terms 
of w defined by 

w 2 = lim lim — — — - (1) 

t^t c - h^o M 2 £ d 

(where \ is the magnetic susceptibility and £ the second-moment correlation length), and 
of the ratios Vj = gj/w^~ 2 . In two dimensions M is known exactly (see e.g. |L7|]) while 
X and £ can be derived using the fact that the two-point function satisfies a Painleve dif- 
ferential equation in the critical limit |[L8|| . Therefore one can easily obtain high-precision 



estimates of w. The ratios Vj are estimated analyzing the LT expansion (to 0(u ) where 
u = e _4//T ) of the free-energy on the square lattice in the presence of an external constant 
magnetic field ||19||. In three dimensions good estimates of w 2 have been obtained from the 



analysis of its LT expansion p0|j20[| and from Monte Carlo simulations [ 21 1 . Field-theoretic 



calculations of w [22,23] are less precise, but perfectly consistent. In order to estimate the 



2 



constants Vj we consider their e-expansion, that can be derived from the e-expansion of the 
equation of state 0-26]. A constrained analysis of these e-series using the corresponding 



two-dimensional results allows us to obtain estimates of v% and V4 with satisfactory accuracy. 
It is straightforward to obtain corresponding estimates for gj = w^ 2 Vj. 

The paper is organized as follows. In Sec. [TJ we introduce our notation and give some 
general formulae for the LT expansion of the effective potential near the coexistence curve. 
In Sec. [ITJ we compute w and Vj for various values of j for the two-dimensional Ising model. 
In Sec. [IV] we present our analysis of the e-expansion of Vj, and we compare our results 
with other approaches. In Sec. [V] we discuss the effective potential for the O(N) model 
in the large- N limit. We find logarithmic terms that have not been previously predicted. 
We discuss some possible interpretations. In App. [A| we report the LT expansion of the 
susceptibilities that are used in the computation of Vj with j < 6 for the two-dimensional 
Ising model. 



II. EXPANSION OF THE EFFECTIVE POTENTIAL NEAR THE 

COEXISTENCE CURVE 

The effective potential is related to the free energy of the model. Indeed, if M = (0) is 
the magnetization and H the magnetic field, one defines 

F{M) = MH-^ log Z(H), (2) 

where Z(H) is the partition function and the dependence on the temperature T is always 
understood in the notation. We will be interested in the behaviour of JF(M) near the 
coexistence curve t = T — T c < 0, H = 0. For the Ising model, if Mq = lim^_».o+ M(H), for 
M > Mq, we can expand J-{M) in powers of M — Mq. Explicitly we can write 

00 1 

F{M) = F(M ) + £ -a 3 (M - M y. (3) 

It is useful to express the coefficients aj in terms of renormalization-group invariant quanti- 
ties. We therefore define a renormalized magnetization 

2 Z(t,H = Q-fM(t,Hf 
* = x( t,J7 = 0) (4) 

where £ is the second-moment correlation length 

2 _ 1 Jdx x 2 G(x) 
5 ~ 2d Jdx G(x) ' 1 ' 

and G(x) is the connected two-point function. In terms of (p we can rewrite 

HV) = HVv) + \m\ V - + E m d -^-^- 9j { V - ^y. (6) 
z j=3 3- 
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Here m = l/£ and gj are functions of t only. In field theory ip is nothing but the zero- 
momentum renormalized field. For t — > 0— the quantities gj approach universal constants 
(that we indicate with the same symbol) that represent the zero-momentum n-point renor- 
malized coupling constants. An even simpler parametrization can be obtained if we introduce 

_ V M 

z = — -TT- {<) 

lf M 



m d 



Then we have 

F(<p)-F(<p ) = ^B(z). (8) 
The function B(z) has the following expansion 

B(z) = hz-iy + Y,U( z -iy, (9) 

j=3 J' 



where 



1 = ^-r do) 



The advantage of the expansion (^|) is that its coefficients Vj are expressed only in terms of 
zero-momentum quantities, while in Eq. (^) the correlation length is also present. 

We mention that a more natural expansion of B(z) that takes into account its parity 
properties would be 



Biz) = - 
v ; 4 



l,o ^ 1 



[Z 



i) 2 + Ett^ 2 -i) j 



2^ n?- 



(11) 



where the coefficients Vj can be easily related to those of the expansion @, i.e. to the 
constants Vj. For example ^3 = (3 — v 3 )/2, etc.... 

Since z oc |t| _/3 M, the equation of state can be written in the form 

\dM) t 11 dz y 1 

where we have used the hyperscaling relation /3(1 + 8) — dv = 0. Eq. (|i~2| ) can be exploited 
to derive B(z) from the equation of state that is written in the form 

H = M 5 f(x) (13) 

with x = tM' 1 ^ . The scaling function f(x) is usually normalized so that /(— 1) = 0, 
/(0) = 1. The function B(z) can be obtained from 



dB{z) 



h z s f (xoz- 1 ^) , (14) 



dz 

where the normalization constants ho and xo are fixed by the requirement that 
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B(z) = ~(z-lf + 0[(z-l) 3 ]. (15) 

Since (3 > and the function f(x) is regular at x = and nonzero, Eq. ([14]) implies that 
B(z) ~ z 5 " 1 " 1 for z — > oo. 

The coefficients Vj of the expansion of B(z) around z = 1 can be related to the LT critical 
limit of combinations of connected Green's functions evaluated at zero momentum 

Xj= E (mtM-'Jfa-iMvj))* (is) 



X2,--.,Xj 



In the LT critical limit 



r> +v 3 , (17) 

Xi M 2 xlM 2 

~)d ^' 

'"^ ^ ( } 

XeM 4 , 1 cM3 M 4 , , v^M 4 inK X4xlM 4 , ln v 4 M 4 



15^ + lO^ 4 — - 105^i — + 105^- ^v 6 , (20) 



-y5 ^.6 y6 -v/7 -y8 

AAA A A 



etc., where M = xi and x = % 2 - 



III. THE TWO-DIMENSIONAL ISING MODEL 

In the critical limit the two-point function of the two-dimensional Ising model, both in 
the symmetric and broken phase, satisfies a Painleve differential equation E8|. Therefore, 
w can be calculated by an appropriate numerical integration, obtaining 

w = 0.72906.... (21) 

The coefficients Vj of the expansion of B(z) are not known. Good estimates of the first 
few Vj can be obtained from the analysis of their LT expansion. The basic reason is that 
the leading correction to scaling is analytic, since the subleading exponent A is expected 
to be larger than one (see e.g. Ref. p7|| and references therein). In particular the available 



exact calculations [18[] for the square-lattice Ising model near criticality have revealed only 
analytic corrections to the leading power law. Therefore the traditional methods of series 
analysis should work well. 



In order to estimate the ratios Vj we used the results of Ref. ||19|| . They report the 
expansion of the free energy in the presence of an external magnetic field to 0(u 23 ) where 
u = e _4//T , for the square-lattice Ising model. Using Eqs. ( |T7l , pOD one can easily obtain 
the corresponding expansion for Vj. In our analysis we used several types of approximants, 
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Pade, Dlog-Pade and first-order integral approximantsQ (for a review on the resummation 
techniques see for example Ref. |28||), constructed from the series of Sj = w~ 2(j ~ 2 ' ] vj that 
have the form Yh=i c i u% - Estimates of Vj are then obtained by evaluating them at u c = (v^2 + 
l) -2 . In App. |A| we report the 0(u 23 ) series of the zero-momentum connected correlation 
functions Xj (with j < 6) we used in our analysis. 
We obtained Q 

v 3 = 33.011(6), (22) 
v A = 48.6(1.2), (23) 
v 5 = 7.69(2) x 10 4 , (24) 
v e w -2.17 x 10 7 , (25) 
v 7 w 9.9 x 10 9 , (26) 
v 8 ss -6.2 x 10 12 . (27) 

Notice how the absolute values of Vj increase rapidly with j. Using the value of w reported 
in Eq. (|2l|) we obtain: 

g 3 = wv 3 = 24.067(4), (28) 

g A = w 2 Vi = 25.8(6). (29) 



In Ref. |T(J the estimates t> 3 = 33.06(10) (f 3 is denoted there by R 3 ) and g 4 f=3 25 (with a 
large apparent uncertainty) were obtained from the analysis of the LT series published in 
Ref. |29| for the square, triangular and honeycomb lattices. These results are fully consistent 
with our analysis. The higher precision we achieved is essentially due to the longer series we 
considered, indeed for the square lattice Ref. [29] reports series to 0(u u ). We also mention 



the estimate g 3 = 23.9(2) reported in Ref. [|1^] that has been obtained by a Monte Carlo 
simulation for £ 



< 



IV. THE THREE-DIMENSIONAL ISING MODEL 

For the three-dimensional Ising model accurate estimates of w can be found in 
the literature. They are obtained using various approaches, from the analysis of low- 



1 Given a rath order series, we considered the following quasi-diagonal approximants: [l/m] Pade 
and Dlog-Pade approximants with I + m > n — 2 and l,m > f — 2; [m/l/k] first-order integral 
approximants with m + l + k + 2 = n and [(n — 2) /3j — 1 < m, I, k < \(n — 2) /3] + 1. 

2 As estimate from each class of approximants, i.e. Pade, Dlog-Pade and first-order integral 
approximants, we considered the average of the values at u = u c of the non-defective approximants 
using all the available terms of the series. As estimate of the error we considered the square root 
of the variance around the estimate of the results from all the non-defective approximants listed 
in footnote |l|. The results from Pade, Dlog-Pade and first-order integral approximants were then 
combined leading to the estimates shown in the text. 
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temperature expansions 
culations 



1Q|J20|| , from Monte Carlo simulations | 2l| , and field-theoretic cal 



22,23,14]. The low-temperature expansion of w 2 can be calculated to 0(u 21 ) on 



the cubic lattice using the series published in Refs. [pO| , pT| . In the analysis of this series we 



used the Roskies transform in order to reduce the systematic effects due to confluent sin- 
gularities. Indeed, since the leading correction-to-scaling exponent A ~ 1/2, the systematic 
error due to the leading non-analytic correction may be relatively large. The analysis using 
Pade, Dlog-Pade and first-order integral approximants (see footnotes |IJ and 0) to the Roskies 
transformed series of w 2 leads to the estimate w 2 = 4.75(4). This is in good agreement with 
the estimate of w 2 that can be obtained from the results of Ref. |TIJ, w 2 = 4.71(5) P3| . The 
result of the Monte Carlo simulation reported in Ref. [21], w 2 = 4.77(3), is consistent with 
these estimates. Therefore, we believe that 



In Refs 



w = 2.18(1) 
|2 



(30) 



a perturbative approach in fixed dimension 
The expansion parameter 



should be a reliable estimate. 

d — 3 is exploited to study the LT phase of the Ising model, 
is u = 3w 2 . Its critical value is determined computing the zero of the corresponding /3- 
function that is known to three loops. The estimate of w 2 one obtains from this expansion is 
consistent with Eq. (|30]), but less precise due to the relatively small number of known term 
in the /^-function. 

In order to calculate the first few v j in three dimensions, we consider their e-expansion. 
The e-expansion of Vj can be derived from the e-expansion of the equation of state that is 
known to 0(e 3 ) for the Ising model p4|-|26[ 



Eqs. ([14] ) and ([IB]) , one obtains 

v 3 = 3 + -e H e 2 + 

3 2 18 

9 43 2 
^4 = 3 + -e + -e- + 



*'5 



Using the procedure described in Sec. |T| cf. 



989 



A C(3)' 



324 




1944 


12 


6 


1601 


A 


C(3) 


648 + 


4 


2 


} 11 X 


+ 


5C(3) X 


4 






12167 




65A 


162 


+ 


2 



0(e% 



6 3 + 0(e 4 ), 



0(e 4 ), 



34C(3) e 3 + 0(e 4 



(31) 
(32) 
(33) 
(34) 



etc., where A = ip'(l/3)/3 - 2n 2 /9 w 1.17195, and C(3) w 1.20206. 

Since the e-expansion is asymptotic, it requires a resummation to get estimates at d = 3, 
i.e. e = 1. Assuming, as usual, its Borel summability, the analysis of the series is performed 
using the method proposed in Ref. pi[ , that is based on the knowledge of the large-order 
behaviour of the series. Given a quantity R with series 



R(e) 



k=0 



we have generated new series R p (a, b; e) according to 

) = ^2B k (a,b)i dtt" 



Rp(a, 



k=0 



u(etf 



[l-u(et)Y 



(35) 



(36) 
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where 



u( X ) = (37) 

Here a = 1/3 is the singularity of the Borel transform. The coefficients Bk(a,b) are deter- 
mined requiring that the expansion in e of R p (a, b; e) coincides with the original series. For 
each a, b and p an estimate of R is simply given by R p (a, fe; e = 1). We follow Refs. p0|,|3| 



in order to derive the estimates and their uncertainty. We determine an integer value of b, 



b opt , such that 



-R3O, b opt ; e = 1) « R 2 (a, b opt ; e = 1) (38) 



for a < 1. & opt is the value of b such that the estimate from the series to order 0(e 3 ) is 
essentially identical to the estimate from the series to order 0(e 2 ). In a somewhat arbitrary 
way we consider as our final estimate the average of R p (a, b;e= 1) with — 1 < a < 1 and 
—2 + b opt <b<2 + b opt . The error we report is the variance of the values of R 3 (a,b;e = 1) 
with — 1 < a < 1 and |_& op t/3 — lj < b < [46 opt /3 + 1]. A discussion of the reliability of 
these error estimates can be found in Refs. pO|J^]. 

When the coefficients of the e-expansion have all the same sign, as it is the case for Vj, 
it is convenient to consider and analyze the series of the inverse. The analysis of the 0(e 3 ) 
series of vj 1 gives 

v 3 = 6.05(27), (39) 
Vi = 17.2(6.6). (40) 

The analysis of the e-series of the coefficients Vj with j > 4 is extremely unstable and does 
not provide reliable estimates. 

These results can be improved performing a constrained analysis that exploits the accu- 
rate two-dimensional estimates of V3 and U4, cf. Eqs. (|22|) and (p3|). The idea is the following. 
Assume that exact or approximate values of the quantity at hand are known for some dimen- 
sions di < 3, with di belonging to the expected analytic domain in d. Then one may use the 
polynomial interpolation between the values d = 4 and d = di as zeroth order approximation 
in three dimensions. If the interpolation is a good approximation one should find that the 
series which gives the deviations has smaller coefficients than the original one. Consequently 
also the errors in the resummation are reduced. The idea to constrain the e-series analysis 



was employed in Ref. |35] to improve the estimates of the critical exponents of the Ising and 



self-avoiding walk models; in Ref. |3(| it was used in the study of the two-point function; in 



Refs. [pO^PI it was successfully applied to the study of the small-renormalized-field expansion 
of the effective potential in the symmetric phase. For quantities defined in the LT of the 
Ising model, the analytic domain in d should contain the value d = 2. Thus we can use 
the rather accurate two-dimensional estimates of v% and V4 to constrain the analysis of the 
e-series. In the study of the effective potential in the symmetric phase also the exact results 
in d — 1, were used, because in that case the analytic domain is expected to extend up to 



d = [|0 
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Assuming that the ratios Vj are analytic and sufficiently smooth in the domain 4 > d > 2 
(that is < e < 2), one may perform a linear interpolation between d = 4 and d = 2, and 
then analyze the series of the difference. For a generic quantity i? one defines 



R(e) 



R{e) - i? ex (e = 2) 



and a new quantity 



# imp (e) = i? ex (e = 2) + (e - 2)R(e) 



(41) 



(42) 



where R ex {£ — 2) is the exact value of R for e = 2. New estimates of i? for e = 1 can 
be obtained applying the resummation procedure we described above to R(e) and then 
computing i? imp (l). 

The constrained analysis of the e-series of vj 1 leads to the following results: 



V3 = 5.99(5 + 0), 
t» 4 = 15.8(1.4 + 0.1) 



(43) 
(44) 



where the second error is obtained varying the two-dimensional estimate within one error 
bar. The new results are in good agreement with the estimates (|39| ) and (fiDD , but have a 
smaller uncertainty. Using the estimate (|30D of w we can compute the three- and four-point 
zero-momentum renormalized coupling constants: 



93 

9a 



13.06(12), 
75(7). 



(45) 
(46) 



The errors of g 3 and g 4 are calculated considering the errors of w, v 3 , and t> 4 as independent. 

Let us compare our results with available estimates obtained using other methods. Table | 
presents a summary of all the available (as far as we know) estimates of t>3, #3 and #4. Field- 
theoretic estimates of V3 are obtained in Refs. ||,|5| using the parametric representation 
of the equation of state (in Table | we refer to this approach by PR). In this approach 
the e-expansion and the d = 3 (^-expansion (i.e. expansion in powers of the HT zero- 
momentum four-point renormalized coupling) are used to estimate the first few coefficients of 
the expansion in powers of 9 of the function h{9) characterizing the parametric representation 
of the equation of state. From h(8) one can derive many universal ratios of quantities 



defined at zero momentum such as Vj. Ref. JTO] presents an analysis of the low-temperature 
expansion on the cubic, b.c.c. and f.c.c. lattices, obtaining t> 3 = 6.47(20), and g 4 ~ 85 
with a large apparent uncertainty. These estimates are in good agreement with our results 
(fJ3|) and ([5|). We also mention Ref. where the effective potential in the broken phase 
is determined from an approximate solution of the exact renormalization-group equations. 
The resulting estimates are w 2 = 5.55 and g% = 15.24. It is very hard to obtain good 



estimates by Monte Carlo simulations. Ref. JT5[ reports some results obtained with a Monte 
Carlo simulation with £ < 8. They obtain g 3 = 16.0(5), and g 4 — 3^3 ~ 750 (the precision 
of the data does not allow us to extract an estimate of g 4 ). Ref. [[LlJ reports estimates of 
#3 and g 4 that are obtained from the probability distribution of the average magnetization, 
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whose numerical data are analyzed by assuming a sixth-order polynomial approximation of 
the effective potential. In Table | we report the data of Ref. ]nj corresponding to the largest 
correlation length £ ~ 7 and lattice 74 3 (taken from the 5th column of Table 2 as suggested 
by the author). While for g 3 there is a substantial agreement, the estimate of g 4 turns out 
to be quite larger than ours. 



V. THE EFFECTIVE POTENTIAL IN THE LARGE- N LIMIT 

The physics of the broken phase of O(N) models with N > 1 is very different from that 
of the Ising model, because of the presence of Goldstone modes. In this case the effective 
potential !F((p) is not analytic for (p — > (p^. General renormalization-group arguments predict 

where c is a constant and p a new exponent that is conjectured to be given exactly by 
p = djid — 2). The conjecture is based on the following argument. The exponent p is 
related to the behaviour of the longitudinal susceptibility \l along the coexistence curve. 
From Eq. (f|7|) it is easy to derive 

Xl = ^~H*-MM (48 ) 

for H — > 0. On the other hand the singularity of xl for H — > is governed by the zero- 
temperature infrared-stable fixed point ||T7j-|59||. This leads to the exact prediction 

XL 

and therefore p = d/(d — 2). The asymptotic behaviour ( f47D has been checked in the large- N 
limit to leading and subleading order |!7| . At leading order Eq. (pTFD holds for all <p and not 
only in the limit (p — > ifo. 

The nature of the corrections to the behaviour ( WJfy is less clear. Setting x = tM^ 1 ^ 
and y = HM~ S , the conjecture is that 1 + x has the form of a double expansion in powers 
of y and y^ d ~ 2 ^ 2 near the coexistence curve |H],[n],[3^] , i.e. for y — ► 



H d ' 2 ~\ (49) 



1 + x = Cl y + C2 y^/ 2 + d lV 2 + d 2 y 2 ~ t/2 + d 3 y 2 -' + ... (50) 

where e = 4 — d. Here we assume, as usual, that x = — 1 corresponds to the coexistence 
curve. With this normalization z = (—x)~^. Thus, for x — > —1, 1+x ~ {z — l)/(3 with 
corrections that are analytic in (z — 1). Therefore, also z — 1 should have an expansion of 
the same form. 

Notice that this expansion predicts that, in three dimensions, z — 1 has an expansion 
that is analytic in powers of y 1 ^ 2 . As a consequence, the effective potential would be analytic 
in z — 1, with -B(^) ~ (z — l) 3 for z —*■ 1. 

These ideas can be verified in the large-A limit using the equation of state to order 
0(1/N) (next-to-leading in 1/N) reported in Ref. P?| . Setting u = 1 + x, we find that, in 
generic dimension d with 2 < d < 4, the function f(x) defined in Eq. (|13|) has an expansion 
of the form 
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f(x) = u 2 '^ 




E E c m ^ 2m/(rf - 2)+n ) + o (- 



-1 m=l 



\N 2 



(51) 



which is consistent with the expansion fl5"0|). Generically, the coefficients CyyiYi 3X6 Sill gular 
for those values of d such that 2m/ (d — 2) is an integer, i.e. for 

2T7X 

2 < d = 2 H < 4, for < m < n, m, n e IN. (52) 

n 

However, in these cases, also the coefficients Cj are singular. A careful analysis shows that 
for these special values of d logarithmic terms should appear. We have verified analytically 
their presence for d = 4 — 2/n, that corresponds to m = n — 1 in Eq. (|52]) . Explicitly in 
three dimensions (n = 2) we find 



f(x) = uS 



1 + 1 (f x (u) + ]oguf 2 (u)) + 0(N~ 2 ) 



(53) 



The functions fi(oo) and f2^) have a regular expansion in powers of u. In particular 



23 1 1 

f 2 (u) = --uj 2 - -(128 + vr V + —(2000 + 17tt V + 0(uj 5 ). (54) 
o lb loO 

Notice the absence of a linear term in /2(a>). This is due to the fact that c\ \ and C\ in Eq. 
(pl|) are not singular for d — > 3. 

The logarithmic terms we find are not predicted by Eq. (|50|) . In particular Eq. ( |53[) is 



incompatible with Eq. (pDD, that, as we already remarked, predicts an analytic expansion 
in powers of uj. The presence of logarithms in the expansion for values of d arbitrarily near 
4 casts some doubts on the asymptotic behavior (|50"D that has been derived essentially from 
an e-expansion analysis. 

What do we learn for the behaviour of the critical effective potential near the coexistence 



curve for finite values of ? A possible interpretation of the results is that the expansion (|53|) 



holds for all values of A^ and thus Eq. (|50D is correct apart from logarithms that are present 
for some special values of d. The reason of their appearance is however unclear. Neverthless, 
it does not necessarily contradict the conjecture that the behavior near the coexistence curve 
is controlled by the zero-temperature infrared- stable Gaussian fixed point. In this case 
logarithms would not be unexpected, as they usually appear in the reduced temperature 
asymptotic expansion around Gaussian fixed points (see e.g. Ref. [Q). Anyway, since 
they appear also for dimensions d > 3, they cannot be related to the presence of marginal 
operators, say (f 2 ) 3 in three dimensions. 

There is another possible interpretation of the result (|53|) that we mention for complete- 
ness. In the large- A" expansion, logarithmic terms are usually the signal of the presence of 
A^-dependent critical exponents. For instance consider a term of the form 

A(N)tu a{N) + B(N)to m , (55) 

with cr(N) = m + a±/N + 0(N~ 2 ) and m integer. For large values of N, expanding A(N) = 
A + Ai/N + 0(N~ 2 ) and analogously B(N), we obtain 
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(A + B )uj m + 1 [{A x + B x )u m + ^u^loga;] + 0(N- 2 ). 

This expansion would reproduce the behaviour (|53"D, provided A + -Bo 
dimension Eq. flSq) would be compatible with the expansion (|5lD if 



(56) 



0. In generic 



2A- 



+ m-2A; + 0(iV 



-In 



(57) 



for some integer k < m/2. Notice that the expansion in powers of 1/N of the coefficients 
A(N) and B(N) would be discontinuous in d, since A(N) and B(N) would be of order 1/N 
in generic dimension and of order 1 whenever <t(oo) is an integer. A similar phenomenon has 
been found for the HT specific heat [43, OH, and in the expansion of the Callan-Symanzik j3- 
function near the critical point [J2U|J . This interpretation of the results is quite natural in the 
framework of the large- N expansion. However it is not clear which exponents should appear. 
In particular, if the singular behaviour is controlled by the zero-temperature infrared-stable 
Gaussian fixed point, we would not expect any iV-dependent exponent. We think that these 
questions deserve further investigation. 
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APPENDIX A: TWO-DIMENSIONAL LT SERIES 



In this appendix we report the LT expansions of the connected correlation functions Xj 
for the the square-lattice Ising model with nearest-neighbor interactions. They are derived 



from the expansion of the free energy in the presence of an external magnetic field [15]. This 
series reproduces the exact results for the free energy and the magnetization in the absence 
of a magnetic field [|5LI7| 




l + w) 4 -16w(l-w) 2 sin 2 # 



1/2 



(Al) 
(A2) 



where u = e _4,/T . This check is particularly effective as all the terms of the expansion 
contribute in the limit H — > 0. 

The expressions for the susceptibilities, cf. Eq. fll6|), up to terms of order 0(u 24 ) are the 
following: 

x = U 2 (A + 32m + 240m 2 + 1664m 3 + 11164m 4 + 73184m 5 + 472064m 6 + 3008032m 7 
+18985364m 8 + 118909888n 9 + 740066448m 10 + 4581660832n n 
+28237063308M 12 + 173353630848m 13 + 1060674765568u 14 
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+6470624695296u 15 + 39370663086596m 16 + 238993166711328m 17 
+1447734754083760m 18 + 8753312020985216m 19 + 52833859249062236u 20 
+318401517346022240m 21 ) (A3) 

X3 = u 2 (-8 - 128m - 1648m 2 - 17216m 3 - 161608m 4 - 1407040m 5 - 11607872m 6 
-91914944m 7 - 704714056m 8 - 5264048832u 9 - 38484337136u 10 
-276313973760m 11 - 1953615567368m 12 - 13630566361728u 13 
-94009039294720m 14 - 641820330787712m 15 - 4342587216599624m 16 
-29147168513086720m 17 - 194228811717309808m 18 - 1285893394538102720u 19 
-8463157853682576776m 20 - 55401729118133053376m 21 ) (A4) 

Xa = m 2 (16 + 512m + 11232m 2 + 174464m 3 + 2268016m 4 + 26007680m 5 
+272557376m 6 + 2666887552m 7 + 24721765328m 8 + 219362198656m 9 
+1877373570336m 10 + 15586471370752U 11 + 126094785804720m 12 
+997562812048896m 13 + 7739640438305152m 14 + 59027832457279488m 15 
+443396861865190928u 16 + 3285761490274027008m 17 + 24053962791473093728m 18 
+174163313822385408896m 19 + 1248496572536156377712m 20 

+8868728488366239392384m 21 ) (A5) 
X5 = u 2 (-32 - 2048m - 76864u 2 - 1763072m 3 - 31561504m 4 - 474188032m 5 
-6284096768m 6 - 75680753408m 7 - 845304372256u 8 - 8883063406848m 9 
-88761722644544m 10 - 850124864318976m 11 - 7853055915373856m 12 
-70313539072037376m 13 - 612653077144892416m 14 - 5211746985723756032m 15 
-43403126238486164768m 16 - 354661581988366266880m 17 
-2849039732654299216960m 18 - 22536604900114543443200m 19 

-175792207497016658107424m 20 - 1353835890235723182430976m 21 ) (A6) 
Xe = m 2 (64 + 8192 u + 531840 m 2 + 17944064 u 3 + 440964544 u 4 + 8650801664 m 5 
+ 144516779264 m 6 + 2136034381312 m 7 + 28672334567744 m 8 
+356014396523008 m 9 + 4144706261854848 m 10 + 45707724350236672 m 11 
+481293469566545088 u 12 + 4869609919772817408 m 13 + 47583052907024977408 u u 
+450914521123508084736 m 15 + 4158352425040048310336 m 16 
+37427697740234745151488 m 17 + 329592281476942345924480 m 18 
+2845676740202773130026496 m 19 + 24132653400815900244380096 m 20 
+201334455541709269472268800 m 21 ) (A7) 

We should note that the expansion of x does not reproduce the expression reported in Ref. 
||19|| . The last two terms are slightly different: 

X(our) - x(Ref. [19]) = 52m 22 + 8 72m 23 (A8) 

The discrepancy is probably due to a misprint, since their general series (in \i = e~ 2H and 



u) correctly reproduces Eqs. (|A1| , |A2|) . 
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TABLES 



TABLE I. We report the available results for V3, 53 and 54. The estimates of g% marked by an 
asterisk have been derived by us using the value of V3 of the corresponding line and the estimate 
© of w. 



e-exp. [this paper] 

e-exp. PR g 
d = 3 -exp. PR § 



ERG |4| 

lt m 



V3 

5.99(5) 

6.07(17) 

6.08(6) 

6.47 

6.47(20) 



MC II 



93 

13.06(12) 
*13.2(4) 
*13. 25(14) 
15.24 
13.9(4) 
13.6(5) 



54 



75(7) 



a 85 
108(7) 
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